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I. INTRODUCTION 



Graphene and carbon nanotubes are excellent materials for nanoelectromechanical sys- 
tems (NEMS). Their large stiffness and low density allow to fabricate high frequency mechan- 
ical resonators sensitive to minute variations of mass, force, and charge 1-25 . In addition to 
their outstanding mechanical properties, these materials owe part of their uniqueness to their 
simplest feature: They constitute the ultimate size limit for one- and two-dimensional (ID 
and 2D) NEMSs. Indeed, single-walled nanotubes are ultra-narrow wires whose diameter 
can be as small as 1 nm; graphene, being a sheet of carbon atoms arranged in a honeycomb 
lattice, is the thinnest membrane imaginable. Figure [T] shows mechanical devices made from 
a carbon nanotube and a graphene sheet. 

Owing to their reduced dimensionality, graphene and carbon nanotubes display unusual 
mechanical behavior. We will discuss two examples in the following. The first one concerns 
the force-displacement response. In the well established case of the cantilever of an atomic 
force microscope probe, its position changes linearly with the force it is subjected to. Such 
a linear force-displacement response is also found in doubly-clamped beams (see Fig. [2h and 




FIG. 1: (a) Scanning electron microscopy image of a nanotube resonator. The nanotube was 
grown by chemical vapor deposition over a prefabricated trench between two W/Pt contacts. The 
nanotube is marked by the arrows and the white dotted lines, (b) Scanning electron microscopy 
image of a suspended single-layer graphene sheet with Au contacts. (Adapted with permission 
from Ref.— .) 
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FIG. 2: (a) Schematic of a suspended beam clamped at both ends and subject to a force. (b,c) 
Expected force-displacement characteristics when the displacement d is much lower than the beam 
thickness t (b) and when d is much larger than t (c). 

b). According to continuum elasticity theory, in the latter case this linearity is expected to 
hold only when the displacement is much smaller than the thickness of the bearn^i. Because 
the diameter of a nanotube and the thickness of a graphene sheet are usually both small 
compared to any displacement, force-displacement measurements in nanotubes and graphene 
are expected to be highly nonlinear (such as in Fig. |2fc). This was indeed shown to be the 
case for graphene^. 

Another striking example of an unusual mechanical property originating from the reduced 
dimensionality of these carbon-based resonators concerns the bending rigidity. This quantity 
characterizes the resistance of an object to bending. Figure shows that the lower region of 
a bent beam contracts while the upper region expands (when the beam bends upwards). At 
the microscopic level, the separation between atoms is larger (smaller) in the upper (lower) 
region. This is energetically not favorable, so the beam tends to return to its initial, straight 
configuration. In this respect, the case of graphene is intriguing (Fig. |3]d): because it is 
only one atom thick, the bending rigidity results solely from the energy cost of changing the 
angle between the p z orbitals of carbon atoms^. Thus the microscopic origin of the bending 
rigidity of graphene differs from the one of standard materials. This book chapter focuses 
on graphene and nanotube NEMS resonators and on the peculiar mechanical properties that 
their reduced dimensionality entails. 
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From a practical point of view, the low dimension and the corresponding ultra-low mass of 
nanotubes and graphene present a great advantage in such experiments as inertial mass sens- 
ing and the exploration of the quantum regime of a macroscopic mechanical degree of free- 
dom. The working principle of mass sensing is simple. Electromechanical resonators can be 
described as harmonic oscillators with an effective mass m (close to the real mass of the res- 
onator volume), a spring constant k, and a mechanical resonant frequency /o = l/2n • \Jk~Jm. 
Mass sensing consists in monitoring the shift in f Q induced by the adsorption of atomic 
species onto the resonator. The reason for the high mass sensitivity of nanotube resonators 
is that the mass of a nanotube is ultra-low, so even a tiny amount of atoms deposited on the 
nanotube makes up a significant fraction of the total mass. Whence, it has been possible 
to achieve a mass sensitivity of about 1 zg = 10 -21 g with a nanotube resonator- - —, which 
surpasses the sensitivity achieved using resonators based on other materials^. Note that 
the frequency shift in inertial mass sensing depends not only on the adsorbed mass but also 
on the position of the adsorbate. By exploiting the nonlinear mechanics of nanoscale res- 
onators, the position can in principle be determined by a measurement at the fundamental 
resonant frequency by using a multi-frequency excitation scheme^. 

Their low mass makes carbon nanotubes and graphene very promising for the study of 
their motion in the quantum regime. This in practice requires preparing the resonator mode 
(oscillator) close to its quantum ground state. What makes such experiments a technological 
feat is that the amplitude of the zero-point motion is typically very low (x zp = h / 4tt f Q m) , 
and therefore difficult to detect. In this context, graphene and nanotubes offer the immediate 
advantage of a very low mass that renders the zero-point motion larger than in heavier 
resonators. A ballpark figure for x zp is typically 1 — 10 fm for resonators microfabricated 
from semiconducting and metallic materials and operated in the 1 — 10 MHz range^i^ 2 -. By 
contrast, x zp is expected to be 1 — 10 pm for graphene resonators with similar resonant 
frequencies. Not only does a large zero-point motion make detecting the motion in the 
quantum limit easier, it also allows for an enhanced coupling to other degrees of freedom, 
such as the photons of a superconducting resonator— ^ and the qubit in a Cooper-pair 
box—. This coupling to nonlinear elements such as two level systems is needed to observe 
quantum dynamics; interestingly, however, graphene resonators are nonlinear at very small 
oscillation amplitudes, opening the possibility to study quantum dynamics without external 
components^. We note that the quantum dynamics of nonlinear vibrations was investigated 
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FIG. 3: (a) Arrangement of the atoms of a beam in a straight (top) and bent (bottom) geometry, 
(b) Arrangement of the p z electronic orbitals of the carbon atoms of a graphene sheet in a straight 
(top) and bent (bottom) geometry. The bending energy is given by the energy cost for the electronic 
orbitals to rotate. 

in the past by means of spectroscopy measurements of the vibrations of impurities in solids 35 . 

II. DETECTION OF THE MECHANICAL VIBRATIONS 

Even though carbon nanotubes and graphene are outstanding materials for NEMS res- 
onators, what first hindered progress in this research field is the difficulty to detect their 
mechanical vibrations. Nanotubes and graphene are so small that it is challenging to apply 
the detection schemes used for larger mechanical resonators 36 . The first measurements were 
carried out by transmission and scanning electron microscopy 1 ' 4 as well as field emission 3 , 
but these detection schemes are suitable neither for sensing nor for experiments in the 
quantum limit. Over the last few years, efforts have been made to develop new detection 
schemes based e.g. on optical interferometry 1 ^, atomic force microscopy^ 1 ^, electrical^, and 
capacitive^ 3 ^ readouts. A technique that has become very popular is the mixing technique. 
It essentially consists of measuring the electrical current flowing through the graphene sheet 
or the nanotube. This detection scheme is very practical since it can be implemented in var- 
ious experimental setups, such as cryostats and ultra-high vacuum chambers. The method 
was first used in experiments on microfabricated resonators^ and was later on adapted to 
nanotube resonators by the McEuen group at Cornell 5 -. The mixing technique was recently 
improved by the group in Lyon= and we will discuss this improved version. 

The resonator is actuated electrostatically by applying an oscillating voltage V FM with 
amplitude V AC to the source electrode, which creates an alternating potential between the 




5 




FIG. 4: (a) Schematic of the device layout and the actuation/detection setup, (b) Schematic of the 
frequency response of the mixing current. The separation between the two minima corresponds to 
the resonance width A/, equal to fo/Q. (Adapted with permission from Ref.— .) 

resonator and the gate electrode (Fig. H^). Central to the technique is the fact that this 
driving voltage is frequency modulated (FM) to take the form V FM = V AC cos[27r/t + 
Ia/Il ■ sin(27r/i)f:)], where / is the carrier frequency and /a is called frequency deviation. 
The motion is detected by measuring the current at frequency fi at the drain electrode, 
-fmix |9Re[xo]/f?/|, where Re[xo] is the real part of the frequency response function of the 
displacement x . The expression |9Re[x ]/9/| can be understood as follows. It contains an 
absolute value b measured with a lock-in amplifier and the output signal is the 

module of the current. The derivative with respect to / is a consequence of the modulation of 
the applied voltage: V FM can be seen as an oscillating voltage whose frequency is modulated 
as fi — f + /a cos(/z,i) (valid over a timescale shorter than Q/fo, where Q is the quality 
factor and /o the resonant frequency of the mechanical resonator). Finally, the origin of 
the real part in |9Re[xo]/9/| is less transparent; as shown in 14 , it comes from the Taylor 
expansion of at small motional amplitude 5x(t) and small drive voltage SV(t), whose 
only component at frequency Jl is contained in the term d 2 I m j X /dxdV ■ 5x5V. 

One advantage of this measurement technique is the frequency conversion. The frequency 
of the measured signal (/&) is about 1 kHz whereas the motion of the resonator can have 
much higher frequencies (typically 10 — 1000 MHz). This frequency conversion is essential 
since the impedance of nanotube and graphene resonators is much larger than 50 Ohm, 
which makes it difficult to measure small high-frequency currents through the resonator. 
A further advantage of this detection scheme is that the mixing current as a function of / 
has a characteristic line-shape (Fig. Hb), which allows us to extract the mechanical quality 
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factor Q in a simple manner: the resonance peak at frequency fo is flanked by two minima 
whose separation is the resonance width Af — fo/Q for a linear harmonic oscillator. Later 
on we will see that the distance between these two minima is also a natural measure of Q 
in a nonlinear oscillator. 

Choosing the right value for the frequency deviation (/a) of the FM technique is crucial for 
a reliable measurement. Namely, one has to ensure that /a is sufficiently small compared 
to the width of the mechanical resonance Af (typically, /a < A//4). Otherwise, the 
measured resonance broadens^^i. (In practice, we measure the dependence of Af on f\; 
the real resonance width is obtained at low /a where Af saturates). 

III. VARIATION OF THE QUALITY FACTOR WITH THE AMPLITUDE OF 
THE MOTION 

Large mechanical resonators, such as those microfabricated from metallic or semicon- 
ducting materials, are usually well described as simple harmonic oscillators. The equation 
of motion is then given by the well-known expression: 

md 2 x/dt 2 = —kx — "fdx/dt + -Fdrive cos(27r ft) (1) 

which interrelates the position x, the velocity dx/dt, and the acceleration d 2 x/dt 2 of the 
oscillator (with effective mass m, spring constant k, and damping rate 7). Here, we show 
that the motion of graphene/nanotube resonators is different from that of larger mechanical 
resonators^. The most striking difference is that the damping depends on the amplitude of 
the motion, in sharp contrast to the behavior of simple harmonic oscillators. 

A surprising experimental fact in nanotube and graphene resonators is that the mechani- 
cal resonance lineshape often broadens as the driving force is increased (-Fd r i ve = C'V g DC V AC 
where C is the derivative of the gate- resonator capacitance with respect to x). In other 
words, the quality factor depends on the amplitude of the motion. This is a novel phe- 
nomenon: indeed, in larger resonators the quality factor is independent of the motional 
amplitude (Q = lixf^raj^). 

Figure [5^ shows an example of the variation of the resonance width (and of the quality 
factor on the right axis) as a function of driving force (which scales linearly with V AC ) in 
a logarithmic plot. The variation of the quality factor is significant and can reach a factor 
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of 50 for some devices. This behavior is robust as it can be observed both in nanotube and 
graphene resonators, both at low temperature (down to 90 mK) and at room temperature, 
and using different detection methods (namely, the FM method described in the previous 
section, the so-called two source technique^, and optical interferometry^) . 

The dependence of the mechanical quality factor on the amplitude of the driving force 
came as a surprise to us. A single mention of this phenomenon appears in^, but its origin 
is not discussed. One possible reason why this effect remained unnoticed before is that the 
variation of the quality factor becomes clearly visible only over a rather large span of driving 
force amplitudes, while the detection of the motion at low drive is challenging. In these 
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FIG. 5: Measurements of the mechanical properties of a resonator made from a nanotube under 
tension. The length of the nanotube is 840 nm and the radius is 2 nm. (a) Resonance width as a 
function of drive V . Squares correspond to 5 K and crosses to 400 mK. The line corresponds 
to equation §5§ with rj = 10 4 kg-m -2 • s -1 (7 = 0). The corresponding quality factor Q is shown 
on the right-hand side scale, (b) Resonance shift as a function of V AC . The line is a solution to 
equation ([3]) with a = 6 • 10 12 kg-m -2 • s -2 . (Adapted with permission from Ref.— .) 
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measurements, care is taken to avoid driving V AC above k^T je in order to prevent electronic 
nonlinear effects or local heating (here &b is the Boltzmann constant, T the temperature, 
and e the electron charge). (On a technical note, the variation of the quality factor does not 
stem from the coupling between electrons and mechanical vibrational^ because the effect 
is not associated to Coulomb blockade and V AC is kept below kftT/e, see Ref.— for a more 
detailed discussion.) 

In order to understand the variation of the quality factor with the drive, it is useful to 
first discuss another behavior. Increasing the driving force also shifts the resonant frequency 
to higher (or sometimes lower) values (Fig. \5jp). This behavior is usually associated to the 
so-called Duffing force -F Duffing = —ax 3 (Ref.—). The latter, together with the spring force, 
can be expressed as 



which indicates that the Duffing term contributes to the restoring force: it makes the res- 
onator stiffer (for a > 0) and increases the resonant frequency. However, the Duffing force 
alone cannot explain the variation of the quality factor with drive, since it does not affect the 
resonance width (in the absence of any bistability), provided that the amplitude of thermal 
fluctuations is smaller than that of the driven motion^. This is readily verified by adding 
the Duffing force to the right hand side of equation (TjQ) and solving for the steady state. The 
resonance width, given by the separation between the two minima in the mixing current as a 
function of / in the FM technique, is indeed independent of the driving force. In addition to 
the Duffing nonlinearity —ax 3 , the other relevant higher-order term in the Newton equation 
for a weakly anharmonic oscillator is the nonlinear damping term —7]x 2 dx/dt (Ref.— ^^). 
We thus get 



A derivation of equation ([3]) based on the Caldeira-Leggett model can be found in4^, in which 
the nonlinear damping force emerges from a nonlinear coupling between the mechanical res- 
onator and a thermal bath of harmonic degrees of freedom (the nonlinearity in the coupling 
results from the anharmonicity of the potential well in which the resonator is confined). 
Other works showed that additional terms of second and third order (x 2 , xdx/dt, (dx/dt) 2 , 
x(dx/dt) 2 , (dx/dt) 3 ) lead in the rotating frame to a renormalization of a and r\ (Ref.— _ — ). 




(2) 




kx — jdx/dt — ax 3 — r]x 2 dx/dt + -Fd r i ve cos(27r ft) 



(3) 
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The term —r]x 2 dx/dt is a damping force since it scales linearly with the velocity. Together 
with —■ydx/dt, it can be expressed as 

- 7 ( 1 + ^-x 2 ^ dx/dt (4) 

The force — rjx 2 dx/dt accounts for a dissipation mechanism that becomes important at large 
motional amplitude. When —'ydx/dt dominates over —t]x 2 dx/dt, which is the case for larger 
mechanical resonators, the resonance width is independent of the driving force and is given 
by Af = ■y/2'iim. In the opposite limits, when the term —jdx/dt can be neglected, 

Af = 0.032m-y/ V /3 ^d/ive (5) 

so that Af oc (y AC ) 2 / 3 . At first sight, the simple relation Q = fo/Af for a linear harmonic 
oscillator is expected to break down in the presence of nonlinearities. However, we show in^ 
that an analogous expression Q = 1.09/o/A/ is recovered in the limit of strong nonlinear 
damping. 

The relation Af oc (y AC ) 2 ^ 3 captures well the experimental data: the line in the double 
logarithmic plot of Fig. |5^i corresponds to a power law with an exponent 2/3. This fit allows 
to extract rj. The agreement between experiment and theory is a strong indication that 
damping is here described by the nonlinear force —r)x 2 dx/dt, instead of the linear force 
—'j dx/dt. 

In the following we discuss two additional experimental facts which further support the 
dominance of the nonlinear damping force in nanotube/graphene resonators. These are 
revealed by studying the hysteretic behavior of the mechanical resonance on one hand and 
parametric excitation on the other. 



IV. HYSTERESIS 



For sufficiently large driving forces, the motional amplitude as a function of the driving 
frequency / develops an asymmetry (dashed line in Fig. [6^, b). This results in bistability 
and hysteresis for certain intervals in / (solid line in Fig. EK, b) 5,21 . The hysteresis is 
intimately related to the resonance shift and also originates from the Duffing force^ 1 ^— . 
An example of a hysteretic resonance lineshape for a graphene resonator is shown in Fig. [6b, 
d. Surprisingly, however, we do not observe a hysteresis in some of our nanotube and 
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FIG. 6: Mechanical properties of a resonator made from a graphene sheet under tensile stress at 
4 K for large driving forces. The length of the sheet is 1.7 fim and the width is 120 nm. (a, b) The 
schematics show the amplitude of motion as a function of driving frequency (solid lines) for both 
sweeping directions. The dashed lines are the solutions of the Duffing equation, (c, d) Frequency 
response of the mixing current at V AC = 500 /xV. The frequency is swept upwards in (c) and 
downwards in (d). (Adapted with permission from Ref.— .) 

graphene resonators. This is particularly intriguing since hysteretic behaviors are always 
observed in larger mechanical resonators. 

The occasional absence of hysteresis is a direct consequence of nonlinear damping. It 
can be predicted from the ratio between the constants of the nonlinear forces a and 77. 
When i]/a > v^3/27r/ , the nonlinear damping is strong enough to keep the broadening 
of the resonance always comparable to or larger than its shift and precludes hysteresis for 
all driving forces^. In the opposite case (when r]/a < y/3/2nf ), a hysteretic behavior is 
expected to emerge. 

In the previous section we saw how r] can be extracted from the drive dependence of the 
quality factor. The parameter a is evaluated by fitting the shift of the resonant frequency as 
a function of V AC to the steady-state solution of equation fl3]) (see solid line in Fig. |5t>). By 
comparing r//a to \/3/2iifo, we can predict the occurrence of the hysteresis. The agreement 
between experiment and theory is good in most cases. This provides the second experimental 
fact that supports the relevance the nonlinear damping force. 
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V. PARAMETRIC EXCITATION 



A whole new range of interesting effects becomes accessible by modulating the parameters 
that enter the Newton equation. For example, nonlinear experiments can be carried out by 
periodically modulating the spring constant k + 5kcos(27if p t). In the case of a pendulum, 
the spring constant can be changed by varying the length of the pendulum's arm (Fig. Wfr)- 
The pendulum can then be brought to resonance by setting f p to be an integer multiple of /o, 
while keeping the driving force to zero. This is what we refer to as "parametric excitation". 

Nanotube resonators are expected to be excellent candidates for parametric excitation 
because k can be modulated with the voltage applied to the gate by a very large amount: it is 
possible to make this modulation larger than in any other mechanical resonators fabricated 
to date (this can be quantified by measuring the gate voltage dependence of the resonant 
frequency, which scales as y/k). Figure [3d shows two clearly resolved resonant modes. Their 
resonant frequency can be tuned with V® c to a large extent. This behavior has been 
attributed to the increase of the elastic tension that builds up in the nanotube as it bends 




V g DC (V) f (MHz) 

FIG. 7: (a) Schematic of a pendulum. By periodically modulating the length Z, one also modulates 
the spring constant k + 5k cos(2ir f p t) with f p an integer number of the resonant frequency /o 
{k = mg/l, where m is the mass of the pendulum and g is the constant of gravitation). In so 
doing, the pendulum can be brought to resonance, (b) Resonant frequency as a function of gate 
voltage (data obtained by measuring Imix versus f and V^^). Two mechanical modes can be seen, 
(c) Parametric excitation signal (dark) obtained by measuring Imix as a function of the detection 
frequency / and Vp. The phase difference between Vp and V^ c (the drive signal) is not kept fixed 
so the measured Imix 

is fluctuating. V^p = 1.4 mV. (Adapted with permission from RefJ^). 
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towards the gate with increasing V® c (Ref.-), in a similar fashion as a bent guitar string 
vibrates at higher frequency. 

In order to realize parametric excitation, we apply an oscillating voltage Vp at a frequency 
2f to the gate electrode (another excitation scheme has recently been proposed by some of 
us^ 7 - and consists in applying a symmetric source-drain voltage). On resonance, this modu- 
lates k at 2/o, thereby achieving parametric pumping of the resonator. The resulting motion 
is detected using the so-called two-source technique 5 , that is, by applying a small voltage 
V s ^° at a slightly detuned frequency (/ — 5f) to the source electrode, and by measuring the 
mixing current 7 m i X at frequency 5f at the drain electrode using a lock-in amplifier. 

Upon increasing the parametric pump excitation above a critical threshold Vpc, the 
nanotube is observed to oscillate without any driving force (Fig. [7b). The resonator becomes 
unstable, so that any fluctuation will activate an oscillating motion that is sustained by the 
parametric drive 4 -. Figure [7b shows mechanical motion in a tongue-shaped region of the 
(/, Vp) space, which is a typical signature of self-oscillation^ 4 ^—. Self-oscillation is observed 
for Vp roughly above Vp t c — 10 mV. In 17 we describe additional measurements (parametric 
amplification measurements) of the same device, which independently yields Vp t c — 10 mV. 

One interesting outcome of this experiment is that the estimation of Vp c allows to extract 
the constant 7 of the linear friction force —'ydx/dt and the associated quality factor Qq = 
27tf m/'-f since below threshold the motional amplitude is so small that nonlinear damping 
can be neglected. Here, we use the simple relation 

Q ° = Vp~c ' dfJdV^ (6) 

where dfo/dV^ is the slope of the V®° dependence of f , and obtain a quality factor 
Qo ~ 1000. This is significantly larger than the quality factor obtained when the parametric 
pump excitation is off (where the quality factor is extracted from the resonance lineshape 
of the resonator driven by -Fd r i ve cos(2tt ft)). There, the quality factor is about 170 — 350 
depending on the driving force (as in section IITT1) . We attribute this difference to the distinct 
damping forces at work: —^dxjdt and —7]x 2 dx/dt. While Vp t c extracted from parametric 
excitation measurements is a measure of 7, the main damping channel is associated to the 
force —t]x 2 dx/dt. This is an additional experimental fact in favor of the nonlinear damping 
scenario in nanotube and graphene resonators. 



13 



VI. DISCUSSION ON THE NONLINEAR DAMPING FORCE 



We have shown that the dynamics of graphene and carbon nanotube mechanical res- 
onators is highly unusual: the quality factor depends on the amplitude of the motion. We 
emphasize however that the quality factor can in principle saturate at low driving force for 
the two following reasons. First, the quality factor is expected to be independent of the driv- 
ing force when the amplitude of thermal vibrations is much larger than that of the driven 
motion, even if the damping is dominated by the —r)x 2 dx/dt force (M. Dykman, private 
communication). Owing to the small mass of nanotubes and graphene sheets, the amplitude 
x t h of the thermal vibrations is large. Using 



yields Xth = 1.4 nm at 300 K for a nanotube with /o = 100 MHz and m = 5 ag (which 
is the mass of a nanotube 1 /im in length and 2.2 nm in diameter). The amplitude of 
the vibrations of driven nanotube resonators was estimated from the measurements to be 
typically 1 — 10 nm at room temperature^ 1 ^. As such, the quality factor is expected 
to rapidly saturate at 300 K at low driving forces (even if the linear therm —'ydx/dt is 
negligible) . 

Naturally, a second mechanism that leads to the saturation of the quality factor at low 
drive is the crossover from a nonlinear to a linear damping regime. Provided that Af and 
the resonance shift are much smaller than f Q , the standard definition of Q is still warranted 
and reads Q = 2nE/AE where AE is the mechanical energy dissipated over one oscillation 
period and E is the corresponding stored energy when the resonator is driven resonantly. 
We obtain 



which shows that the quality factor depends on the driving force at large Xq whereas it 
becomes a constant at low xq (here, xq is the maximum steady state amplitude on resonance). 

Figure M shows that the quality factor of a nanotube resonator measured at 300 K satu- 
rates at low driving force. The dashed line is a fit to the model that assumes the occurrence 
of a crossover from the nonlinear to the linear damping regime using 7 and 77 as fit parame- 
ters (we note that the device is different from that in Fig. [5]). We emphasize that the data 
can alternatively be captured by a model that assumes that the saturation is due to the 




(7) 




(8) 
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FIG. 8: Measurements of the mechanical properties of a resonator made from a nanotube at 300 K. 
Resonance width and quality factor as a function of driving force measured with the FM technique 
(filled squares) and 2-source technique (hollow squares), respectively. The solid line represents a 
comparison to equation ([5]) with negligible linear damping (77 = 2.5 • 10 3 kg-m~ 2 • s -1 , 7 = 0), the 
dashed line is obtained with a finite 7 (77 = 1 • 10 3 kg-m~ 2 • s _1 , 7 = 1.9 • 10~ 14 kg-s -1 ). (Adapted 
with permission from Ref.— .) 

thermal vibrations. More work is required to identify the proper contributions of each of 
these two effects. 

The theoretical foundations of damping rest on Newton's Principia. Damping has been 
successfully described by the linear damping force —^dx/dt for all mechanical resonators 
studied so far in vacuum. Remarkably, this picture holds for resonators whose dimensions 
span many orders of magnitude down to a few tens of nanometers. Reducing dimensions to 
the atomic scale using graphene and nanotube resonators, our work demonstrates that the 
simple linear damping scenario ceases to be valid and damping is often better described by 
incorporating —i]x 2 dx/dt. 

This finding has profound consequences in light of the fact that several predictions re- 
garding quantum and sensing experiments assume that damping is linear. Such exper- 
iments, if carried out with nanotube or graphene resonators, would indeed have to be 
interpreted within the framework of nonlinear damping. These include studies of the 
quantum-to-classical transition^, the cooling efficiency^, the mass resolution 55 , and the 
force sensitivity 56 . 

Our results not only provide new insight into the dynamics of nanotube and graphene 
resonators, they also hold promise for technological applications. Our control over the res- 
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FIG. 9: Frequency response of the mixing current for a graphene resonator. Each panel corresponds 
to a given gate voltage. The solid lines are fits with Q = 10 5 . Measurements are carried out at 
90 mK with V AC = 8 fxV. The integration time of the lock-in amplifier is 300 ms. /a = 1 kHz 
and similar results are obtained with /a = 500 Hz. The resonance width is typically 1.5 kHz (and 
occasionally 300 Hz smaller or larger). (Adapted with permission from Ref.— .) 

onance width allows us to improve the mechanical quality factor. In order to achieve larger 
Q-factors, we simply lower the driving force until the motion becomes barely detectable. 
For this, it is convenient to select the value of V® c for which the detection signal is largest. 
In so doing, we measured a quality factor of 10 5 for a graphene resonator at 90 mK (see 
Fig- E]), which is the largest Q ever reported in a graphene resonator—. Larger quality fac- 
tors enable better force sensing. We obtained a force sensitivity of 2.5 aN/vTIz using a 
nanotube resonator operating at 100 mK. This is within a factor of five of the best sensi- 
tivities reported using microfabricated resonators operating at their ultimate limit set by 
thermal vibrations^- 1 ^. 
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The microscopic origin of the nonlinear damping is still to be determined, but it could 
be related to such diverse phenomena as phonon tunneling into the supports^^, sliding 
at the contacts, nonlinearities involving phonon-phonon interactions, or contamination in 
combination with geometrical nonlinearities^. It would be of considerable interest to ex- 
perimentally study the dependence of the nonlinear damping force on contamination, the 
clamping configuration, and the suspended length. Theoretical work on the microscopic 
nature of nonlinear damping should also prove valuable^ 1 ^. 

VII. COUPLING MECHANICAL VIBRATIONS TO ELECTRON TRANSPORT 
IN A NANOTUBE RESONATOR 

Lastly, we show that strong electromechanical nonlinearities arise when a nanotube res- 
onator operated in the Coulomb blockade regime vibrates in concert with single electron 
hopping. 

Figure [10] shows the mechanical properties of a nanotube resonator at 4 K. Both the 
resonant frequency and the quality factor oscillate as a function of the voltage applied to 
the backgate electrode. These measurements show that the dynamics of nanotubes can be 
widely tuned by external electric means {V g DC ), which is very advantageous for applications. 

To understand the origin of the oscillations of the resonant frequency and of the quality 
factor, it is useful to first consider the electrical properties of the nanotube. The conductance 
oscillates with V g DC in a way that is typical of the Coulomb-blockade regime and with the 
same V g DC period as the oscillations of f and Q (see Ref.— ). This correlation indicates that 
the mechanical motion is affected by Coulomb blockade. 

Charge transport through nanotubes in the Coulomb-blockade regime is a well-studied 
phenomenon. Upon sweeping V g DC , the conductance oscillates as the charge gdot residing 
on the nanotube increases stepwise (Fig. [TT]) . Furthermore, in our experiment the nanotube 
behaves as a mechanical resonator. The capacitance C g between the nanotube and the gate 
oscillates in time as 5C g = C' g x (with C' g the derivative of C g with respect to displacement 
x). As a result, the nanotube is periodically charging and discharging by the amount <5gdot = 
SC g V g DC (as shown in Fig. [TH . 

The oscillating charge on the nanotube results in a shift of the resonant frequency of the 
resonator. When electrons tunnel onto and out of the nanotube, the relevant electrostatic 
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FIG. 10: (a) Resonant frequency and (b) quality factor as a function of gate voltage V® c for a 
nanotube resonator at 4 K. The nanotube length is 1 /im and the diameter is ~ 1.1 nm. (Adapted 
with permission from Ref.— .) 



force acting on the nanotube can be expressed as^ 

F P -- 



C ,yDC 

(gdot - q c ) 



a 



dot 



(9) 



where q c = —C g V® c , called control charge, is introduced for commodity (this force comes 
from F e = 0.5C^(Vdot — V® ) 2 with the electrostatic potential of the nanotube dot Vdot = 
(<7dot ~ ?c)/Cdot)- Because of the repeated charging and discharging of the nanotube by the 
amount <5gdot, equation flU]) results in a spring force, F e = —5k ■ x. This force modifies the 
spring constant of the nanotube resonator by 5k. As a result, the resonator softens and the 
resonant frequency gets reduced by the amount 6 fo = (fo/2) ■ (5k /k). 

Besides, charging and discharging the nanotube causes damping of the mechanical motion. 
Indeed, the charge <5gdot has to flow across the tunnel barrier at the nanotube-electrode 
interface, and the dissipated energy of the charge (after the tunnel process) is supplied from 
the mechanical resonator. 

This damping and the resonant frequency are expected to oscillate with V® c , since the 
amplitude of the £gdot oscillations depends on V® c . The amplitude is large (low) when gdot 
increases sharply (weakly) with V g DC ; see Fig. CCD Our model captures the data well as 
discussed in^. We note that the background signal of the resonant frequency in Fig. [TUh 
decreases continuously over the full V® c sweep: this is attributed to the increase of the 
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V e (=-^ v 9 /e ) 

FIG. 11: Schematic of the charge on the nanotube dot as a function of the control charge (q c = 
—C g V® c ) in the Coulomb-blockade regime. Vibrations cause the charge on the nanotube to 
oscillate (with amplitude 5qdot and A^dot for low and large motional amplitude, respectively). 
(Adapted with permission from Ref.— .) 

elastic tension that builds up in the nanotube as it bends towards the gate upon increasing 
V® c —. Second, the (oscillating) dips of the resonant frequency in Fig.fTUa are not symmetric, 
which may be caused by the stepwise change of the electrostatic potential when sweeping 
V g DC (Ref.i2). 

As we have seen, the coupling between mechanical vibrations and charge transport results 
in an electrostatic force acting on the nanotube^ 

-^electro ^electro lelectrodx / dt (10) 

where the constants k e i ectro and 7 e i ec tro can be tuned with V® c , which is very practical for fu- 
ture use. 7 e iectro is so large that the associated damping force dominates over the —r\x l dxjdt 
force. Equation (Q is valid for low Sq^ot- The situation changes for large mechanical oscilla- 
tion amplitude. The charging and discharging of the dot becomes highly non-trivial during 
one oscillation cycle (see Ag dot in Fig. [TT|) . As a result, -F electro is expected to depend non- 
linearly on x and dx/dt (in contrast to the linear dependences at low oscillation amplitude) 
and .Feiectro has to be solved numerically^. 

We found evidences for nonlinearities in F e i ec tro- Figures fT2h.b show two resonance peaks 
at 1.5 K (mixing current as a function of driving frequency measured using the two-source 
technique^). The resonance in Fig. [T2h has a surprising lineshape -it splits into two peaks, 
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FIG. 12: Mixing current I m ix as a function of driving frequency for two different gate voltages 
at 1.5 K. A second, narrow peak appears in (a). The inset shows the nanotube conductance as a 
function of gate voltage V g ■ AV® C is measured from the maximum of the Coulomb-blockade 
peak. The nanotube is the same as in Fig. [TUJ (Adapted with permission from RefJ^.) 

one broad and the other narrow. The narrow peak appears only at gate voltages which yield a 
Coulomb blockade peak in the conductance (inset of Fig. fl~2b). Strikingly, its width narrows 
down as the driving force increases. We show that the double-peak structure arises from 
the interplay between two nonlinear phenomena, one associated with the nonlinear force 
-^electro and the other with the nonlinear detection mechanism. Indeed, the nonlinearity in 
the detection mechanism due to Coulomb blockade can affect the lineshape of the frequency- 
dependent amplitude x(f) in such a way that the lineshape of I m ix(f) develops a double 
peak structure (the reason for this is that for sufficiently large motional amplitude I m i X 
decreases when x increases)^. If the x(f) lineshape were a Lorentzian, the double-peak 
structure in I mix would be symmetric. However, the nonlinear force -F e i ec tro renders the x(f) 
lineshape strongly asymmetric, which results in two peaks in I m ix(f), one broad and the 
other narrow-^. 

Overall, these results show that the coupling between mechanical motion and electron 
transport in nanotube resonators can be made so strong that the associated force acting 
on the nanotube becomes highly nonlinear in its displacement and velocity. This strong 
coupling originates from the reduced dimension of nanotubes. Larger and heavier mechanical 
resonators are evidently much less sensitive to the motion of individual tunneling electrons^. 

We note that this electro-mechanical coupling also influences charge transport in the 
Coulomb blockade regime, a topic that has attracted much interest^— . Another active field 
of research consists in harnessing this coupling to cool the mechanical motion of resonators 
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to the phononic ground stated— , as well as controling the dynamics of large cantilevers of 
scanning probe microscopes^— . 



VIII. CONCLUSION 



Graphene and carbon nanotubes are very interesting systems for studying resonant me- 
chanical behavior. They constitute the ultimate size limit of one and two-dimensional 
NEMSs: nanotubes are wires with a diameter of about 1 nm and graphene is a membrane 
that is only one atom thick. Because of their reduced dimensionality, graphene and carbon 
nanotubes display unusual mechanical phenomena, among which strong nonlinearities are 
ubiquitous. In this chapter, we reviewed several types of nonlinear behavior. We first dis- 
cussed an unprecedented scenario where damping is described by a nonlinear force. This 
scenario is supported by several experimental facts: (i) the quality factor varies with the 
amplitude of the motion as a power law whose exponent coincides with the value predicted 
by the nonlinear damping model, (ii) hysteretic behavior (of the motional amplitude as a 
function of driving frequency) is absent in some of our resonators even for large driving 
forces, as expected when nonlinear damping forces are large, and (iii) when we quantify 
the linear damping force (by performing parametric excitation measurements) we find that 
it is significantly smaller than the nonlinear damping force. We illustrated this chapter 
with measurements on nanotube resonators but we also observed the experimental facts 
(i) and (ii) on graphene resonators (see Ref.— ). We then reviewed parametric excitation 
measurements on nanotube resonators, an alternative actuation method which is based on 
nonlinear dynamics. Finally, we discussed experiments where the mechanical motion is cou- 
pled to electron transport through a nanotube. The coupling can be made so strong that 
the associated force acting on the nanotube becomes highly nonlinear with displacement 
and velocity. Overall, graphene and nanotube resonators hold promise for future studies on 
classical and quantum nonlinear dynamics. 
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